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Preface 


The  sediment  impact  assessment  for  Manitou  Springs,  CO,  was  conducted 
at  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  at  the 
request  of  the  U.S.  Army  Engineer  District,  Albuquerque  (SWA).  It  was 
prepared  in  support  of  the  Manitou  Springs  Reconnaissance  Report. 

This  investigation  was  conducted  during  the  period  September  to 
November  1994  in  the  Hydraulics  Laboratory  of  WES,  under  the  direction  of 
Messrs.  Frank  A.  Herrmann,  Jr.,  Director  of  the  Hydraulics  Laboratory; 
Richard  A.  Sager,  Assistant  Director  of  the  Hydraulics  Laboratory;  Dr.  Larry 
L.  Daggett,  Chief  of  the  Waterways  Division  (WD),  Hydraulics  Laboratory; 
and  Mr.  Michael  J.  Trawle,  Chief  of  the  Math  Modeling  Branch  (MMB), 

WD.  The  project  engineer  was  Dr.  Ronald  R.  Copeland,  MMB.  Authors  of 
the  report  were  Dr.  Copeland  and  Ms.  Lisa  C.  Hubbard,  MMB. 

During  the  course  of  this  study,  close  working  contact  was  maintained  with 
Mr.  Bruce  Beach,  SWA,  who  served  as  coordinating  engineer,  providing 
required  data,  technical  assistance,  and  review 

At  the  time  of  publication  of  this  report,  Director  of  WES  was 
Dr.  Robert  W.  Whalin.  Commander  was  COL  Bruce  K.  Howard,  EN. 


The  contents  of  this  report  are  not  to  be  used  for  advertising,  publication , 
or  promotional  purposes .  Citation  of  trade  names  does  not  constitute  an 
official  endorsement  or  approval  of  the  use  of  such  commercial  products. 


Conversion  Factors, 
Non-SI  to  SI  Units  of 
Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  units 
as  follows: 


Multiply 

By 

To  Obtain 

acre-feet 

4,046.873 

cubic  meters 

cubic  feet 

0.02831685 

cubic  meters 

cubic  yards 

0.7645549 

cubic  meters 

inches 

2.54 

centimeters 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  kelvins1 

feet 

0.3048 

meters 

miles  (U.S.  statute) 

1.609347 

kilometers 

pounds  (force)  per  square  foot 

47.88 

Pascals 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  meter 

square  miles 

2.589998 

square  kilometers 

To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings,  use  the  following 
formula:  C  =  (5/9)(F  -  32).  To  obtain  Kelvin  (K)  readings,  use:  K  =  (5/9)(F  -  32)  +  273.15. 

1  Introduction 


Location 

The  city  of  Manitou  Springs  is  located  in  central  Colorado,  just  west  of 
Colorado  Springs,  in  a  narrow,  steep-walled  mountain  canyon  that  contains 
Fountain  Creek.  The  town  was  developed  as  a  resort  around  the  natural 
mineral  springs  in  the  area,  and  remains  a  popular  tourist  destination.  Over 
the  years,  residential  and  commercial  development  has  encroached  into  the 
natural  streambed  of  Fountain  Creek.  The  creek  is  further  constricted  by 
several  bridges  and  culvert  crossings  and  in  some  locations  by  major  roads 
paralleling  the  creek.  These  alterations  to  the  natural  creek  channel  have 
increased  the  risk  of  flooding  and  channel  instability  in  Manitou  Springs.  The 
study  area  extends  from  the  downstream,  or  eastern,  city  limit  of  Manitou 
Springs,  near  Columbia  Road,  upstream  to  a  location  about  2,000  ft1 
upstream  from  Park  Avenue.  A  location  and  vicinity  map  is  shown  in 
Figure  1. 


Watershed  Characteristics 


The  drainage  area  of  Fountain  Creek  at  the  eastern  city  limit  of  Manitou 
Springs  is  about  103  square  miles.  The  headwaters  are  located  in  the  moun¬ 
tains  of  the  Rampart  Range.  The  drainage  basin  includes  Pikes  Peak,  which 
rises  to  an  elevation  of  14, 1 10  ft.2  The  elevation  in  Fountain  Creek  at  the 
downstream  study  limit  is  about  6,060  ft.  The  watershed  is  typically  moun¬ 
tain  terrain  with  steep  rough  slopes.  Fountain  Creek  is  a  perennial  stream, 
fed  from  glacial  snowpacks  and  springs.  A  watershed  map  is  shown  in 
Figure  2. 

Vegetation  varies  throughout  the  drainage  basin.  In  areas  above  the 
timberline,  vegetation  is  sparse.  In  the  timber  and  woodland  areas  aspen,  oak 
brush,  spruce,  juniper,  and  native  grasses  are  found. 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  found  on 
page  vi. 

2  All  elevations  (el)  cited  herein  are  in  feet  referred  to  the  National  Geodetic  Vertical  Datum 
(NGVD). 
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Figure  1.  Location  and  vicinity  map 


Watershed 

Boundary 


Four  tributaries  join  Fountain  Creek  in  the  study  reach,  as  listed  in  the 
following  tabulation: 


Creek 

Drainage  Area 
square  miles 

Fountain  Creek  upstream  from  Manitou  Springs 

70 

Ruxton  Creek 

17.8 

Williams  Canyon  Creek 

2.7 

Sutherland  Creek 

5.2 

Black  Canyon  Creek 

2.6 

Fountain  Creek  downstream  from  Manitou  Springs 

103 

Ruxton  Creek  flows  in  a  northeasterly  direction  to  join  Fountain  Creek  at 
Manitou  Springs.  Like  that  of  Fountain  Creek,  the  watershed  is  typically 
steep  mountain  terrain  with  steep  rough  slopes.  Elevations  range  between 
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14,110  ft  at  Pikes  Peak  to  6,330  ft  at  the  confluence  with  Fountain  Creek. 

Two  water  supply  reservoirs  are  in  the  watershed:  Lake  Moraine  and  Big 
Tooth  Reservoir.  The  last  half  mile  of  the  channel  is  constricted  by  bridges 
and  culverts  as  it  traverses  residential  and  commercial  developments.  The  last 
400  ft  of  Ruxton  Creek  upstream  from  its  confluence  with  Fountain  Creek  are 
completely  covered. 

The  Williams  Canyon  Creek  watershed  lies  north  of  Manitou  Springs. 
Williams  Canyon  Creek  flows  in  a  southerly  direction  from  its  headwaters  in 
the  Rampart  Range  to  join  Fountain  Creek  just  below  the  confluence  of 
Ruxton  Creek.  Elevations  in  the  watershed  range  from  9,400  ft  to  6,320  ft. 
The  original  creek  channel  has  been  replaced  by  underground  conduits 
through  the  city  of  Manitou  Springs  from  the  mouth  of  Williams  Canyon  for 
about  1,000  ft  to  Fountain  Creek. 

Sutherland  Creek  flows  in  a  northeasterly  direction  and  joins  Fountain 
Creek  just  upstream  from  the  U.S.  Highway  24  crossing.  The  Sutherland 
Creek  watershed  lies  south  and  east  of  Ruxton  Creek  watershed.  Elevations 
range  from  10,800  ft  to  6,220  ft.  The  confluence  of  Sutherland  Creek  with 
Fountain  Creek  has  been  relocated  due  to  highway  construction.  In  its  present 
location  the  creek  is  severely  constricted  upstream  from  its  confluence  with 
Fountain  Creek  by  a  small  culvert  that  passes  under  Manitou  Avenue. 

Black  Canyon  Creek  lies  directly  east  of  Williams  Canyon  Creek  and  flows 
in  a  southerly  direction  to  join  Fountain  Creek  just  downstream  from  the 
U.S.  Highway  24  crossing.  Elevations  in  the  watershed  range  from  8,200  to 
6,210  ft. 


Valley  and  Streambed  Characteristics 


Fountain  Creek’s  valley  through  Manitou  Springs  upstream  from  U.S. 
Highway  24  has  an  average  width  of  about  400  ft,  and  the  channel  varies  in 
width  between  12  and  50  ft.  Through  town  the  channel  is  severely  constricted 
by  bridges  and  building  foundations.  Several  channel  stabilizers  with  1-  or  2- 
ft  drops  have  been  constructed  across  the  streambed.  Downstream  from  U.S. 
Highway  24,  the  valley  width  averages  about  600  ft,  and  the  channel  width 
varies  between  50  and  100  ft. 

The  composition  of  the  streambed  varies.  In  some  locations  the  bed  is 
covered  with  cobbles  and  boulders,  in  other  locations  sand  and  gravel.  In 
some  locations  there  are  abundant  quantities  of  rubble  on  the  bed.  Although 
not  observed  during  the  field  investigation,  an  earlier  report  (U.S.  Army 
Engineer  District  (USAED),  Albuquerque,  1974)  stated  that  bedrock  outcrops 
are  also  present. 
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Climate  and  Precipitation 


In  the  mountainous  areas,  precipitation  varies  widely  over  relatively  short 
distances,  and  much  of  the  total  precipitation  at  these  higher  elevations  is  in 
the  form  of  snow.  In  the  Fountain  Creek  watershed,  snowmelt  seldom  pro¬ 
duces  floods  except  when  augmented  by  rainfall.  Average  annual  precipitation 
varies  between  24  in.  at  Lake  Moraine  and  13  in.  at  Colorado  Springs.  Most 
of  the  flood-producing  storms  occur  between  May  and  August.  The  physical 
features  of  the  watershed  are  all  conducive  to  a  rapid  concentration  of  runoff 
resulting  in  flash  floods  characterized  by  high  peak  flows,  moderate  volumes, 
and  short  durations. 

Temperatures  vary  widely  because  of  altitude  differences.  Mean  annual 
maximum  temperatures  vary  between  47  and  63  °F.  Mean  annual  minimum 
temperatures  vary  between  24  and  35  °F.  Recorded  extreme  temperatures 
vary  between  -37  and  100  °F. 


Purpose  of  the  Study 

The  purpose  of  the  sediment  impact  assessment  is  to  identify  the  magnitude 
of  sediment  problems  that  might  be  associated  with  proposed  flood-control 
projects  for  Manitou  Springs  and  to  recommend  appropriate  sedimentation 
studies  for  the  next  level  of  planning  study.  For  a  reconnaissance  level  plan¬ 
ning  study,  this  may  be  accomplished  by  evaluating  channel  stability  using  a 
sediment  budget  approach.  Bed  material  sediment  yields  for  specific  channel 
reaches,  and  for  project  and  existing  conditions  are  compared.  When  there 
are  significant  differences  in  calculated  sediment  yields  between  one  reach  and 
another  or  between  existing  and  project  conditions,  a  sedimentation  problem 
has  been  identified.  The  sediment  impact  assessment  also  provides  an  inven¬ 
tory  of  available  sediment  data  and  recommends  data  collection  programs  if 
appropriate. 

The  sediment  impact  assessment  also  includes  estimates  of  sediment  yield 
for  purposes  of  sizing  reservoirs  and/or  debris  basins.  These  estimates  are 
approximate  in  the  reconnaissance  level  study  due  to  the  lack  of  detailed 
watershed  data.  Data  collection  requirements  for  more  detailed  sediment  yield 
studies  are  identified  in  the  sediment  impact  assessment. 

At  the  time  of  this  study  no  specific  flood-control  plans  had  been  formu¬ 
lated,  so  comparison  of  existing  and  project  plans  was  limited.  The  sediment 
impacts  of  a  plan  to  divert  all  flows  greater  than  500  cfs  just  downstream 
from  Park  Avenue  and  then  return  the  diverted  flows  downstream  from 
U.S.  Highway  24  was  evaluated. 
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Approach 


Sediment  transport  and  sediment  yields  for  the  sediment  budget  approach 
were  determined  using  techniques  based  on  uniform  flow  assumptions.  This 
technique  is  recognized  to  be  untenable  for  final  design,  especially  in  streams 
such  as  Fountain  Creek  with  widely  variable  characteristics.  However,  it  is 
deemed  adequate  to  identify  the  magnitude  of  sediment  problems  and  to  eval¬ 
uate  relative  impacts  for  a  variety  of  proposed  flood-control  project  plans. 
Hydraulic  parameters  for  the  sediment  budget  were  determined  from  reach- 
averaged  values  taken  from  the  HEC-2  backwater  model  (U.S.  Army 
Engineer  Hydrologic  Engineering  Center  (USAEHEC)  1990).  The  SAM 
hydraulic  design  package  (Thomas  et  al.,  in  preparation)  was  used  to  calculate 
sediment  transport  and  the  bed  material  sediment  yield. 

There  is  no  generally  accepted  method  to  calculate  total  sediment  yield. 

In  order  to  obtain  a  range  of  probable  sediment  yields  for  a  1  percent  chance 
exceedance  hydrograph  and  for  average  annual  conditions,  sediment  yields 
were  calculated  using  methods  developed  for  watersheds  with  somewhat 
similar  conditions.  This  approach  provides  a  range  of  volume  requirements 
for  construction  of  reservoirs  or  debris  basins.  Further  study  would  be 
required  if  the  sediment  storage  option  is  pursued  in  the  next  level  of  study. 
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2  Channel  Stability 


Channel  Geometry 

The  Albuquerque  District  developed  a  HEC-2  (USAEHEC  1990)  back¬ 
water  model  using  cross  sections  taken  from  topographic  mapping  with  a  scale 
of  1  in.  =  400  ft  with  5-ft  contour  intervals,  modified  somewhat  using  field 
surveys.  The  model  extended  from  the  eastern  city  limits  of  Manitou  Springs 
near  Columbia  Road  for  about  14,200  ft  upstream  to  the  western  city  limits. 
This  model  was  used  in  the  sediment  impact  assessment  study  to  obtain  reach- 
averaged  values  for  important  hydraulic  variables. 

Plots  of  calculated  average  boundary  shear  stress,  channel  discharge, 
channel  velocity,  and  Froude  number  for  the  1  percent  and  50  percent  chance 
peak  exceedance  discharges  are  shown  in  Figures  3-10.  Distances  in  plots  are 
measured  in  feet  upstream  from  Columbia  Road.  These  critical  hydraulic 
parameters  vary  significantly  through  the  study  reach.  This  is  due  partially  to 
poor  channel  definition  at  some  of  the  cross  sections,  but  is  due  primarily  to 
numerous  constrictions  by  culverts,  bridges,  and  channel  encroachment  by 
structures.  This  variability  suggests  a  wide  range  in  streambed  characteristics 
and  channel  stability  conditions  in  the  study  reach,  and  makes  it  difficult  to 
apply  simple  reach-averaged  analysis  techniques  to  Fountain  Creek.  An 
adequate  sediment  analysis  would  require  a  HEC-6  numerical  simulation  that 
accounts  for  nonuniform  water  and  sediment  discharges  (USAEHEC  1993). 

The  reach-averaged  hydraulic  variables  of  width,  depth,  slope,  and  velocity 
must  be  considered  approximate,  because  the  topographic  mapping  scale  used 
to  define  the  cross  sections  is  not  detailed  enough  to  obtain  sufficiently  accu¬ 
rate  channel  definition.  In  addition,  the  backwater  model  did  not  include  the 
effect  of  channel  stabilizers,  many  of  which  cause  a  drop  in  the  channel  invert 
of  between  1  and  2  ft.  This  effect  may  be  significant  at  lower  discharges. 
Future  sedimentation  studies  should  be  conducted  using  new  topographic 
mapping  at  a  scale  of  1  in.  =  100  ft,  which  will  provide  much  better  geo¬ 
metric  definition.  However,  cross  sections  developed  from  the  new  topo¬ 
graphic  mapping  will  still  require  refinements  based  on  field  surveys. 
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3.  Boundary  shear  stress,  1  percent  chance  peak  exceedence  discharge 


Figure  4.  Boundary  shear  stress,  50  percent  chance  peak  exceedence  discharge 
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Figure  5.  Channel  discharge,  1  percent  chance  peak  exceedence  discharge 


Figure  6.  Channel  discharge,  50  percent  chance  peak  exceedence  discharge 
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Figure  9.  Froude  number,  1  percent  chance  peak  exceedence  discharge 


Hydraulic  Parameters  for  Reaches 


The  study  reach  was  divided  into  five  reaches  for  the  sediment  budget 
analysis.  Average  hydraulic  parameters  were  determined  from  the  HEC-2 
output  for  each  reach  using  the  SAM  hydraulic  design  package.  The  SAM 
program  weights  calculated  hydraulic  parameters  by  reach  length.  An  attempt 
was  made  in  the  designation  of  reach  boundaries  to  avoid  cross  sections 
through  bridges  and  immediately  upstream  from  bridges.  The  upstream  reach 
extended  from  Park  Avenue  upstream  for  1,900  ft.  This  reach  was  designated 
as  the  supply  reach.  The  next  downstream  reach,  Reach  2,  extended  from  the 
confluence  of  Ruxton  Creek  downstream  to  Lafayette  Avenue,  and  included 
the  most  commercially  developed  portion  of  Fountain  Creek.  Reach  3  began 
at  Old  Mans  Trail  and  extended  downstream  to  Buena  Vista  Avenue.  The 
channel  in  this  reach  was  also  constricted  at  many  locations.  Downstream 
between  Buena  Vista  Avenue  and  Beckers  Lane  was  Reach  4.  This  reach  was 
relatively  unconstricted.  The  most  downstream  reach  was  Reach  5  between 
Beckers  Lane  and  Manitou  Avenue.  This  reach  was  also  relatively  uncon¬ 
stricted,  but  had  a  milder  average  slope  than  Reach  4.  Profiles  of  Fountain 
Creek  with  reach  designations  are  shown  in  Figures  11-13. 


Hydrology 

The  peak  discharges  and  hydrographs  used  in  this  sediment  impact  assess¬ 
ment  were  taken  from  a  hydrology  report  prepared  by  Albuquerque  District 
(USAED,  Albuquerque,  1994).  Peak  discharges  were  obtained  from  a 
regional  frequency  analysis.  Reported  peak  discharges  are  listed  in  the 
tabulation  on  page  16. 

The  Albuquerque  District  used  the  HEC-1  Flood  Hydrograph  Package 
(USAEHEC  1981)  to  obtain  flood  hydrographs  for  the  Fountain  Creek  water¬ 
shed.  The  model  was  calibrated  to  reproduce  the  peak  discharge  determined 
from  the  regional  frequency  analysis.  The  1  percent  chance  exceedance 
hydrographs  for  three  locations  on  Fountain  Creek  are  shown  in  Figure  14. 

Rainfall  for  various  percent  chance  exceedance  frequencies  for  the  6-hour 
and  1-hour  storms  were  obtained  from  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)  Atlas  2,  Volume  III,  for  Colorado  (Miller,  Frederick, 
and  Tracey  1973).  Depth-area  curves  obtained  from  NOAA  Atlas  2  were 
used  to  determine  the  average  depth  for  each  drainage  basin. 

A  flow  duration  curve  was  developed  to  analyze  average  annual  sediment 
yield.  The  flow  duration  curve  was  based  on  data  from  the  U.S.  Geological 
Survey  (USGS)  Fountain  Creek  near  Colorado  Springs  gage,  which  has 
operated  since  1958.  This  gage  is  located  about  0.5  mile  downstream  of  the 
U.S.  Highway  24  bridge,  and  monitors  a  drainage  area  of  103  square  miles. 
The  flow  duration  curve  developed  from  33  years  of  mean  daily  flow  records 
was  modified  to  account  for  peak  flows.  The  USGS  publishes  all  peak  flows 
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Figure  1 1 .  Water-surface  profiles  for  supply  reach  and  Reach  2 


REACH  3 
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Figure  12.  Water-surface  profiles  for  Reach  3 


Drainage 

Area 

Percent  Chance  That  Discharge,  cfs.  Will  Be  Exceeded  in 

Any  Year 

Percent  Chance 

Location 

miles 

50 

10 

4 

2 

1 

0.2 

Fountain  Creek 
upstream  of 
Ruxton  Creek 

70 

340 

1,400 

3,400 

5,300 

7,300 

13,400 

Fountain  Creek 
downstream  of 
Ruxton  Creek 

88 

340 

1,700 

4,000 

6,500 

9,300 

16,400 

Fountain  Creek 
downstream  of 
Sutherland 

Creek 

103 

340 

1,960 

4,600 

7,000 

10,400 

17,600 

Ruxton  Creek  at 
Fountain  Creek 

17.8 

_1 

370 

940 

1,700 

2,800 

5,600 

Sutherland 

Creek  at 

Fountain  Creek 

5.2 

_  1 

110 

370 

690 

1,090 

2,100 

Williams 

Canyon  Creek 
at  Fountain 

Creek 

2.7 

61 

250 

420 

660 

1,160 

Black  Canyon 
Creek  at 

Fountain  Creek 

2.6 

_1 

60 

220 

410 

650 

1,250 

1  Not  available. 

over  250  cfs  for  this  gage.  Hydrographs  for  these  peak  flows  were  assumed 
to  have  the  same  shape  as  the  1  percent  chance  exceedance  hydrograph  com¬ 
puted  for  Fountain  Creek.  Exceedance  times  for  these  peak  hydrographs  were 
included  in  the  flow  duration  curve  shown  in  Figure  15. 

Albuquerque  District  (USAED,  Albuquerque,  1994)  determined  that  the 
Fountain  Creek  near  Colorado  Springs  gage  has  experienced  peak  flows  that 
are  unrealistically  low  when  compared  to  other  stream  gauges  in  the  region. 
This  suggests  that  the  average  annual  sediment  yields  calculated  with  a  flow 
duration  curve  developed  using  data  from  this  gage  may  produce  values  that 
are  also  too  low. 


Bed  Material  Gradation 


Fountain  Creek  through  Manitou  Springs  has  variable  bed  characteristics. 
In  general,  the  bed  is  composed  of  cobbles  and  boulders  overlying  a  finer  sub¬ 
strate  that  contains  sediment  sizes  from  medium  sand  through  boulders.  A 
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Figure  14.  One  percent  chance  exceedance  hydrograph  for  Fountain  Creek 

variety  of  construction  activities  within  the  channel  have  provided  effective 
control  points.  These  include  concrete  inverts,  log  stabilizers,  and  dumped 
rubble.  These  local  hydraulic  controls  typically  reduce  velocities  upstream, 
resulting  in  a  sand  and  gravel  bed.  In  other  locations  the  bed  is  effectively 
armored  with  larger  gravels,  cobbles,  and  boulders. 

During  a  field  reconnaissance  in  September  1994,  samples  of  the  coarse 
surface  layer  and  of  the  finer  sand  and  gravel  bed  material  were  obtained. 

The  coarse  surface  layer  was  determined  in  Fountain  Creek  about  300  ft 
downstream  from  Park  Avenue.  The  surface  layer  gradation  was  determined 
using  the  Wolman  (1954)  method.  One  hundred  grains  were  randomly  col¬ 
lected  from  the  bed  surface  and  measured.  Fifty  feet  upstream  from  Mayfair 
Avenue  a  bed  sample  was  taken  from  a  point  bar  just  above  the  water  surface 
level  about  midway  along  the  bar.  The  sample  was  collected  with  a  shovel 
about  4  in.  deep  into  the  bar.  The  bed  gradations  are  plotted  in  Figure  16. 
For  purposes  of  this  study,  the  surface  layer  gradation  was  used  to  determine 
critical  shear  stress,  and  the  point  bar  gradation  was  assumed  to  represent  a 
subsurface  gradation  and  was  used  to  calculate  sediment  transport  once  the 
critical  shear  stress  was  exceeded. 
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Figure  15.  Flow  duration  curve,  water  years  1958-1992 

Critical  Shear  Stress 

The  coarse  surface  layer  on  the  bed  of  much  of  Fountain  Creek  serves  to 
prevent  the  underlying  sand  and  gravel  from  becoming  entrained  into  the  flow 
and  thus  protects  the  bed  from  erosion.  However,  at  high  flows  this  armor 
layer  may  be  destroyed,  exposing  the  underlying  material  to  entrainment  and 
transport.  When  this  happens,  sediment  transport  potential  is  greatly 
increased.  To  calculate  a  sediment  transport  rating  curve  for  a  range  of  dis¬ 
charges,  it  is  necessary  to  estimate  the  discharge  at  which  the  armor  layer  is 
no  longer  effective.  Gessler  (1971)  determined  from  flume  studies  that  the 
movement  of  coarse-grained  material  from  the  surface  layer  was  a  probabilis¬ 
tic  process  and  that  a  portion  of  the  coarse  surface  layer  would  remain  effec¬ 
tive  until  the  applied  bed  shear  stress  was  about  twice  the  critical  shear  stress. 
The  critical  shear  stress  can  be  calculated  from  Equation  1 : 


Tc  -  0.047  (ys  y)  d^Q 


(1) 
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Figure  16.  Bed  gradations— surface  and  subsurface 
where 

tc  =  critical  shear  stress 
Ys  =  specific  weight  of  sediment 
7  =  specific  weight  of  water 
d50  =  median  grain  size 

The  surface  layer  sample  from  Fountain  Creek  downstream  from  Park  Avenue 
indicated  an  armor  layer  median  grain  size  of  80  mm  (0.2625  ft).  Assuming  a 
specific  weight  of  sediment  of  165  lbf/cu  ft  yields  a  critical  shear  stress  of 
1.27  lbf/sq  ft. 

Average  boundary  shear  stress  in  the  channel  for  each  cross  section  in  the 
HEC-2  backwater  model  was  calculated  for  a  range  of  discharges  and  com¬ 
pared  to  the  critical  shear  stress  and  twice  the  critical  shear  stress.  These 
comparisons  for  100  and  340  cfs  (50  percent  chance  peak  exceedance  fre¬ 
quency)  are  shown  in  Figures  17  and  18.  At  100  cfs,  boundary  shear  stress  is 
below  the  critical  shear  stress  for  the  majority  of  the  cross  sections  and  below 
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twice  the  critical  shear  stress  for  all  of  the  cross  sections.  Therefore,  at  this 
discharge,  the  armor  layer  was  assumed  to  be  intact  and  the  armor  layer 
gradation  was  used  in  the  sediment  transport  calculations.  At  340  cfs  the 
boundary  shear  stress  in  the  channel  was  greater  than  the  critical  shear  stress 
for  a  majority  of  the  cross  sections,  and  above  twice  the  critical  shear  stress 
for  some  of  the  cross  sections.  Therefore,  at  this  and  greater  discharges  the 
armor  layer  was  assumed  to  be  ineffective  and  the  finer  bed  gradation  was 
used  to  calculate  sediment  transport.  This  analysis  assumes  that  the  bed 
gradation  taken  from  the  point  bar  upstream  from  Mayfair  Avenue  is 
representative  of  the  typical  subsurface  gradation  in  Fountain  Creek 
throughout  the  study  reach.  This  is  a  significant  assumption  that  should  be 
verified  in  a  more  detailed  sediment  sampling  program. 


Sediment  Transport 

Sediment  transport  in  Fountain  Creek  can  be  divided  into  two  categories: 
wash  load  and  bed  material  load.  Wash  load  is  the  finer  sediment  load  that  is 
supplied  by  the  watershed  and  does  not  depend  on  the  hydraulic  characteristics 
in  the  channel.  Wash  load  does  not  exchange  with  the  bed,  although  some 
wash  load  grain  sizes  may  appear  in  the  bed  due  to  entrapment  behind 
localized  obstructions  or  deposits  due  to  receding  flows.  Einstein  (1950)  sug¬ 
gested  that  the  lower  10  percent  of  a  bed  gradation  curve  should  be  discounted 
to  account  for  wash  load  in  the  bed.  Bed  material  load  is  the  load  that 
actively  exchanges  with  the  bed  and  can  be  calculated  given  the  bed  material 
gradation  and  the  hydraulic  characteristics  of  the  channel.  Sediment  transport 
equations  predict  only  the  bed  material  load.  Bed  material  load  is  the  sedi¬ 
ment  load  that  affects  channel  change  and  is  the  important  load  when 
considering  channel  stability. 

There  is  no  generally  acceptable  sediment  transport  equation  for  streams 
characteristic  of  Fountain  Creek.  Sediment  transport  equations  are  notorious 
for  predicting  widely  variable  results.  This  is  especially  true  for  streams  with 
coarse  surface  layers.  The  surface  layer  will  provide  a  variable  degree  of 
protection  to  the  underlying  material  depending  on  the  hydraulic  conditions. 
Upstream  supply  to  any  given  reach  is  also  critical  to  the  development  of  the 
surface  layer  and  the  sediment  transport  rate.  In  Fountain  Creek,  sediment 
transport  rates  will  also  be  highly  variable  due  to  the  variability  in 
cross-section  shape  and  to  the  flashy  nature  of  the  flood  hydrographs.  The 
approach  taken  here  was  to  calculate  sediment  yield  using  several  transport 
functions  and  then  choose  three  that  produced  a  reasonable  range  of  results, 
recognizing  that  quantitative  results  cannot  be  accurately  obtained  using  a 
uniform  flow  approach. 

Much  better  quantitative  results  can  be  obtained  using  the  HEC-6 
numerical  sedimentation  model  (USAEHEC  1993).  This  model  will  account 
for  variable  armoring  and  for  the  variation  in  channel  cross-sectional  areas. 
Careful  adjustment  of  the  model  will  produce  a  more  reliable  estimate  of  bed 
material  transport  through  Manitou  Springs. 
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Bed  Material  Sediment  Yield 


Bed  material  sediment  yield  was  calculated  using  the  flow  duration- 
sediment  discharge-rating  curve  method  described  in  Engineer  Manual  (EM) 
1110-2-4000  (Headquarters,  U.S.  Army  Corps  of  Engineers,  1989).  The 
SAM  hydraulic  design  package  was  used  to  make  the  calculations.  Flow 
duration  curves  were  developed  for  the  mean  annual  flow,  the  50  percent 
chance  exceedance  hydrograph,  and  the  1  percent  chance  exceedance  hydro¬ 
graph.  Sediment  discharge  rating  curves,  required  for  this  method,  were 
developed  for  each  reach  using  three  sediment  transport  functions. 

The  three  sediment  transport  functions  used  in  this  study  were  chosen  from 
eleven  different  sediment  transport  functions  available  in  the  SAM  hydraulic 
design  package.  The  eleven  equations  chosen  for  initial  consideration  have 
been  used  successfully  on  other  sand  and  gravel  bed  streams.  These  equations 
are  divided  into  two  categories:  those  developed  to  calculate  the  bed-load 
component  of  the  bed  material  load,  and  those  developed  to  calculate  the  total 
bed  material  load.  In  order  to  compare  the  sediment  transport  equations,  sedi¬ 
ment  yield  in  the  supply  reach  was  calculated  for  the  1  percent  exceedance 
hydrograph,  the  50  percent  exceedance  hydrograph,  and  for  average  annual 
conditions  using  the  flow  duration  curve.  Results  are  presented  in  Table  A1 . 
There  was  a  wide  range  in  calculated  yields  depending  on  which  equation  was 
used.  Since  there  are  no  available  measured  sediment  concentration  data  to 
circumstantiate  the  sediment  transport  equations,  the  choice  of  equations  was 
based  on  judgment.  Bed-load  equations  were  excluded  because  suspended 
load  may  be  a  significant  contributor  to  the  sediment  load  at  high  flows.  The 
Laursen-Copeland  and  Englund-Hansen  equations  were  excluded  because  cal¬ 
culated  concentrations  at  low  flows  were  unreasonably  high.  A  reasonable 
range  of  calculated  sediment  yields  were  obtained  using  the  Yang,  Ackers- 
White,  and  Toffaleti-Schoklitsch  equations. 

The  following  sediment  yields  were  calculated  from  the  study  reach  using 
the  three  sediment  transport  functions: 


Sediment  Transport  Equation 

Bed  Material  Load,  cu  yd 

Average 

Annual 

50  Percent 

Chance 

Exceedance 

Hydrograph 

. 

1  Percent 

Chance 

Exceedance 

Hydrograph 

Ackers-White 

90 

310 

4,500 

Yang 

150 

660 

10,200 

Toffaleti-Schoklitsch 

1  380 

920 

7,100 

The  calculations  are  based  on  reach-averaged  hydraulic  parameters.  Bed 
material  gradations  were  based  on  the  surface  or  subsurface  sample  as  dis¬ 
cussed  in  the  section,  “Bed  Material  Gradation.” 
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Sediment  Budget 


The  calculated  sediment  yields  for  each  reach  were  compared  to  the  calcu¬ 
lated  sediment  yield  for  the  reach  immediately  upstream  to  determine  if  there 
was  a  trend  for  aggradation  or  degradation.  This  treatment  means  that  the 
most  upstream  reach  is  taken  to  be  the  supply  reach  and  is  assumed  to  be  in 
equilibrium.  Different  sediment  transport  equations  produced  different  calcu¬ 
lated  sediment  yield  quantities  and  therefore  different  amounts  of  degradation 
or  aggradation.  However,  each  equation  predicted  relatively  consistent  values 
in  terms  of  the  percent  of  the  upstream  sediment  load  that  could  be 
transported  through  the  reach.  The  sediment  budget  does  not  account  for  the 
change  in  sediment  transport  capacity  that  would  occur  with  significant  deposi¬ 
tion  or  scour,  and  therefore  the  calculated  aggradation  or  degradation  would 
be  overpredicted.  Calculated  percentages  averaged  from  the  three  sediment 
transport  equations  are  listed  in  the  following  tabulation: 


Reach 

Percent  of  Sediment  from  Upstream  Reach  Passing  this  Reach 

Average  Annual 

50  percent  Chance 
Exceedance 

Hydrograph 

1  percent  Chance 
Exceedance 

Hydrograph 

2 

48 

62 

64 

3 

103 

98 

120 

4 

186 

170 

272 

5 

68 

81 

50 

The  sediment  budget  indicates  that  in  Reach  2,  downstream  from  Ruxton 
Creek  to  Lafayette  Avenue,  between  52  and  36  percent  of  the  incoming  bed 
material  sediment  load  would  be  expected  to  deposit  in  the  channel  depending 
on  the  hydrograph.  This  deposition  is  attributed  to  channel  constrictions  and 
the  resulting  backwater  in  this  reach,  which  passes  through  the  most  commer¬ 
cially  developed  portion  of  Fountain  Creek.  Reach  3,  the  next  reach  down¬ 
stream,  between  Old  Mans  Trail  and  Buena  Vista  Avenue,  can  be  expected  to 
experience  degradation  during  the  1  percent  chance  exceedance  hydrograph, 
but  should  be  stable  in  terms  of  aggradation  and  degradation  during  the 
50  percent  chance  exceedance  hydrograph  and  on  an  average  annual  basis. 
Reach  4,  between  Buena  Vista  Avenue  and  Beckers  Lane,  is  steeper  than 
Reach  3  upstream,  and  has  fewer  constrictions.  As  a  result,  sediment 
transport  potential  is  significantly  greater  in  this  reach,  and  the  sediment 
budget  analysis  indicates  that  the  channel  should  degrade,  especially  during 
the  1  percent  chance  exceedance  hydrograph.  The  most  downstream  reach, 
between  Beckers  Lane  and  Manitou  Avenue,  also  lacks  constrictions,  but  has  a 
milder  slope  than  the  upstream  reach  so  that  aggradation  is  predicted. 
Calculated  sediment  yields  are  presented  in  Table  A2. 
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The  sediment  budget  analysis  indicates  significant  variability  in  sediment 
transport  capacity  in  Fountain  Creek  through  Manitou  Springs.  This  suggests 
channel  stability  problems  even  with  the  existing  channel  conditions.  To 
quantify  the  magnitude  of  the  channel  stability  problem,  better  definition  of 
the  channel  bed  material  and  channel  cross  section  will  be  required.  A  more 
detailed  sediment  study  should  include  a  thorough  channel  bed  inventory, 
identifying  both  surface  and  subsurface  gradations  and  bedrock  outcrops, 
along  Fountain  Creek.  It  is  expected  that  the  bed  material  gradations  will 
vary  longitudinally  through  the  study  reach.  Special  care  should  be  taken  to 
identify  bed  conditions  at  locations  where  degradation  is  predicted.  Locations 
and  reliability  of  channel  bed  stabilizers  should  also  be  determined. 


Impacts  of  Diversion  at  Ruxton  Creek 


One  flood-control  alternative  that  is  being  considered  is  the  diversion  of 
flood  flows  in  excess  of  channel  capacity  downstream  from  Park  Avenue. 

The  diverted  flows  would  then  be  returned  to  Fountain  Creek  downstream  of 
U.S.  Highway  24.  Details  for  this  alternative  had  not  been  developed  in  time 
for  inclusion  in  this  report;  however,  the  impacts  of  the  plan  on  sediment 
transport  can  be  evaluated  using  the  sediment  budget  approach. 

The  diversion  structure  should  be  constructed  on  the  outside  of  a  channel 
bend  so  that  a  minimum  of  bed  material  load  will  be  diverted.  The  design 
should  include  a  side  overflow  weir  that  will  exclude  normal  low  flows  and 
bed  load.  The  exclusion  of  bed  load  from  the  covered  diversion  is  important 
because  deposits  in  the  diversion  conduit  would  create  conveyance  and  main¬ 
tenance  problems.  For  purposes  of  the  sediment  impact  assessment,  it  was 
assumed  that  the  diversion  would  carry  all  flood  flows  greater  than  500  cfs. 

It  was  further  assumed  that  all  the  bed  material  load  would  remain  in  Fountain 
Creek.  Hydrographs  were  modified,  and  the  sediment  budget  calculations 
provided  the  percent  passage  of  bed  material  yield  from  the  upstream  reach 
shown  in  the  following  tabulation. 

The  sediment  budget  analysis  indicates  that  the  diversion  project  will  result 
in  more  deposition  in  Reach  2  downstream  from  the  proposed  diversion  for 
the  1  percent  chance  exceedance  hydrograph  and  for  the  average  annual  basis. 
For  the  50  percent  chance  exceedance  hydrograph,  there  is  no  difference 
because  no  flow  is  diverted.  Significantly  more  degradation  will  occur  in 
Reach  4,  downstream  from  the  diversion  return  point,  for  the  1  percent 
chance  exceedance  hydrograph  and  the  average  annual  basis.  This  is  due  to 
the  reduced  sediment  supply  from  upstream  with  the  implementation  of  the 
diversion  plan.  Increased  maintenance  costs  in  the  existing  channel  should  be 
included  in  the  design  plan  for  the  diversion  structure. 
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Percentage  of  Sediment  from  Upstream  Reach  Passing  this  Reach 
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Approach 

Calculations  of  total  sediment  yield  were  made  for  each  watershed  contri¬ 
buting  sediment  to  Fountain  Creek  within  the  study  reach.  These  calculations 
were  made  to  obtain  estimates  of  sediment  volumes  for  sizing  reservoirs 
and/or  debris  basins.  Estimates  of  total  sediment  yields  were  made  for  both 
the  50  percent  chance  exceedance  flood  and  the  1  percent  chance  exceedance 
flood  and  for  average  annual  conditions. 

The  total  sediment  yield  is  composed  of  both  wash  load  and  bed  material 
load.  As  defined  previously,  the  wash  load  is  the  fine  sediment  that  remains 
in  suspension  without  exchanging  with  the  bed  once  it  reaches  a  channel. 

Wash  load  sources  are  the  watershed  surface,  gullies,  and  the  channel  bed  and 
banks.  The  bed  material  load  is  the  sediment  load  that  actively  exchanges 
with  the  channel  bed  as  it  is  transported  downstream  as  either  suspended  load 
or  bed  load.  The  bed  material  load  capacity  is  determined  by  the  composition 
of  the  bed  and  the  hydraulic  properties  of  the  channel. 

There  is  no  generally  accepted  method  for  calculating  total  sediment  yield. 
Available  techniques  require  measured  sediment  deposition  or  transport  data 
for  adjustments  to  establish  coefficients.  Because  many  factors  affect  the  sedi¬ 
ment  yield,  it  is  generally  necessary  to  have  a  significant  sediment  database  to 
refine  a  technique  to  the  point  where  it  can  be  used  to  make  reliable  predic¬ 
tions.  This  database  does  not  exist  in  the  Manitou  Springs  study  area.  The 
approach  taken  in  this  study  was  to  apply  techniques  that  have  been  used  in 
somewhat  similar  watersheds,  using  limited  available  data,  and  then  draw 
some  general  conclusions  about  the  magnitude  and  uncertainty  of  the  sediment 
yield. 

Sediment  yields  from  steep  mountainous  watersheds  during  a  storm  event 
can  be  estimated  using  the  Los  Angeles  District  Method  (USAED,  Los 
Angeles.  1992).  This  empirical  method  was  developed  using  data  from  water¬ 
sheds  in  the  Southern  California  Coastal  Range.  Sediment  yields  predicted  by 
this  method  represent  sediment  trapped  in  debris  basins  and  consist  primarily 
of  coarser  sediment  sizes.  The  average  annual  sediment  yield  can  be 
predicted  using  the  Pacific  Southwest  Inter-agency  Committee  (PSIAC) 


Chapter  3  Total  Sediment  Yield 


Method  (PSIAC  1968).  This  is  also  an  empirical  method  based  on 
depositional  data  from  the  Pacific  southwest. 


Los  Angeles  District  Method 

The  Los  Angeles  District  Method  is  based  on  a  statistical  analysis  of  mea¬ 
sured  deposition  in  debris  basins,  hydrologic  data,  and  watershed  characteris¬ 
tics.  The  database  for  these  equations  includes  that  of  Tatum  (1963)  plus 
additional  data  collected  from  debris  basins  located  in  the  Southern  California 
Coastal  Range  subsequent  to  Tatum’s  work.  This  method  is  intended  to  esti¬ 
mate  the  debris  yield  from  coastal-draining,  mountainous  Southern  California 
watersheds.  Outside  of  the  recommended  application  area,  careful  adjustment 
of  the  calculated  yields  is  required. 

A  total  of  350  observations  from  80  watersheds  were  used  to  develop  the 
following  regression  equation  used  to  calculate  unit  debris  yield  for  drainage 
areas  up  to  3  square  miles: 


log  Dy  =  0.65  (log  F)  +  0.62  (log  RR)  +  0.18  (log  A) 
+  0.12  (FF) 


(2) 


where 

Dy  =  unit  debris  yield,  cu  yd/square  mile 

P  =  maximum  1-hour  precipitation,  in  inches,  taken  to  two  decimal 
places  and  multiplied  by  100 

RR  -  relief  ratio,  ft/mile 

A  =  drainage  area,  acres 

FF  =  fire  factor 

The  50  percent  and  1  percent  chance  exceedance  regional  1-hour  point 
rainfalls  were  obtained  from  Miller,  Frederick,  and  Tracey  (1973).  Technical 
Memorandum  NWS  HYDRO-40  (NOAA  1984)  was  used  to  reduce  the  point 
precipitation  for  a  given  drainage  area,  which  is  the  maximum  1-hour 
precipitation  required  by  Equation  2.  It  should  be  noted  that  one  of  the 
limitations  of  the  Los  Angeles  District  Method  is  that  the  flood  should  have  at 
least  a  20  percent  chance  exceedance  frequency.  Using  the  50  percent  chance 
exceedance  flood  is  an  extrapolation  of  the  method,  and  results  must  be  used 
with  caution.  The  relief  ratio  is  determined  as  follows: 
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HP  -  LP 


(3) 


RR  = 


L 


Is 


where 

HP  =  highest  point  in  the  watershed  taken  at  the  extension  of  the  longest 
stream,  ft 

LP  =  lowest  point,  in  this  case  the  debris  collection  site,  ft 

Lls  =  length  of  the  longest  stream  up  to  the  highest  point  in  the 
watershed,  miles 

Several  of  the  drainage  basins  considered  in  this  study  have  dams.  Any 
sediment  produced  upstream  of  the  dams  would  not  reach  the  debris  collection 
site.  In  this  study,  the  noncontributing  areas  upstream  from  dams  have  been 
neglected  when  determining  subbasin  drainage  areas.  Fountain  Creek 
upstream  from  Ruxton  Creek  has  several  dams  that  reduce  its  drainage  area 
from  70  square  miles  to  50  square  miles,  and  the  drainage  area  of  Ruxton 
Creek  at  Fountain  Creek  is  reduced  from  18  square  miles  to  12  square  miles. 
The  nondimensional  fire  factor  is  found  from  a  set  of  fire  factor  curves. 

These  curves  are  dependent  on  the  size  of  the  drainage  area.  The  first  curve 
(Figure  19)  is  applicable  for  a  drainage  area  between  0.1  and  3  square  miles. 
The  remaining  curves  (Figure  20)  are  for  a  range  of  drainage  areas  from  3  to 
200  square  miles.  When  the  destruction  of  the  watershed  vegetation  by  fire  is 
not  considered,  the  fire  factor  is  set  to  3,  an  unburned  value.  It  was  con¬ 
sidered  unlikely  that  a  full  or  partial  burn  would  be  followed  by  a  1  percent 
chance  exceedance  flood.  In  this  instance  an  unburned  fire  factor  was  con¬ 
sidered  reasonable.  To  predict  sediment  yield  after  a  total  watershed  burn,  the 
50  percent  chance  exceedance  flood  was  considered. 

The  calculated  debris  yield  from  Equation  2,  in  statistical  terms,  is  the 
expected  value.  Uncertainty  associated  with  the  calculated  result  can  be  mea¬ 
sured  using  the  standard  deviation  of  the  estimate  of  the  expected  value.  The 
standard  deviation  for  Equation  2  is  0.465  (log  Dy).  It  can  be  stated  with 
67  percent  confidence  that  the  “true”  value  of  debris  yield  is  within  one 
standard  deviation  of  the  expected  value.  It  can  also  be  stated  with  95  percent 
confidence  that  the  true  debris  yield  will  fall  within  two  standard  deviations  of 
the  expected  value.  These  statistics  are  based  on  the  data  used  to  develop  the 
regression  equation  and  assume  that  any  calculated  value  comes  from  a  water¬ 
shed  with  similar  geomorphic  and  hydrologic  conditions. 

For  drainage  areas  between  3  and  10  square  miles,  the  following  regression 
equation  was  developed: 
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log  Dy  =  0.85  (log  0  +  0.53  (log  RR)  (4) 

+  0.04 (log  A)  +  0.22 (FF) 


where  Q  is  the  unit  peak  discharge,  cfs/square  mile.  The  unit  peak  discharge 
is  found  by  dividing  the  peak  discharge  by  the  drainage  area. 

The  equation  for  drainage  areas  between  10  and  25  square  miles  is: 

log  Dy  =  0.88  (log  0  +  0.48 (log  RR)  +  0.06 (log  A) 

+  0.20  (FF)  1  ; 


When  the  drainage  area  is  between  25  and  50  square  miles  the  equation  is: 

log  Dy  =  0.94 (log  0  +  0.32 (log  RR)  +  0.14(log  A) 

+  0A7(FF)  ^  ; 
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YEARS  SINCE  100%  WILDFIRE 


Figure  20.  Fire  factor  curves  for  watersheds  from  3.0  to  200  square  miles 

For  larger  drainage  areas  between  50  and  200  square  miles  the  equation  is: 

log  Dy  =  1.02 (log  Q)  +  0.23 (log  RR)  +  0.16(log  A) 

+  0. 13  (FF)  (7) 


A  total  of  187  observations  from  seven  watersheds  were  used  in  the 
development  of  Equations  4  through  7.  The  standard  deviation  was  deter¬ 
mined  to  be  0.242  log  D 

Fountain  Creek  above  Ruxton  Creek  has  a  drainage  area  of  50.3  square 
miles.  Because  this  size  is  very  close  to  the  boundary  determining  the  use  of 
either  Equation  6  or  7,  in  this  study  both  equations  were  used  to  calculate  the 
sediment  yield.  The  difference  in  sediment  yield  obtained  from  using  Equa¬ 
tions  6  and  7  indicates  a  discontinuity  in  the  system  of  equations  and  provides 
another  indication  of  the  uncertainty  associated  with  the  technique. 

In  the  development  of  the  Los  Angeles  District  Method  regression  equa¬ 
tions,  watersheds  with  high  unit  yields  were  used.  Applying  this  technique  to 
areas  with  less  erosion  will  overestimate  debris  yields.  This  led  to  the 
development  of  the  adjustment  and  transposition  (A-T)  factor.  This  factor 
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takes  into  account  the  unquantifiable  geomorphic  and  geologic  parameters  that 
affect  debris  production.  Watersheds  in  the  San  Gabriel  Mountains  of 
southern  California,  which  provided  the  data  for  the  regression  equations, 
would  have  an  A-T  factor  of  1 .  Areas  that  have  less  debris  yield  potential 
would  have  values  less  than  1  and  those  with  a  greater  potential  would  use  a 
value  greater  than  1 .  The  Los  Angeles  District  Method  provides  several 
techniques  for  determining  the  A-T  factor.  The  preferred  techniques  require 
data  from  the  subject  or  nearby  watersheds.  The  required  data  include  mea¬ 
sured  deposition  in  debris  basins  from  storms  with  known  runoff  or  rainfall, 
average  annual  rainfall,  and  sediment  yield.  When  no  debris  yield  data  of  any 
kind  are  available,  a  technique  is  used  that  requires  a  detailed  field  analysis 
identifying  geomorphic  characteristics  of  the  watershed.  As  no  debris  yield 
data  were  available  for  the  area  under  consideration  or  a  neighboring  area,  the 
latter  method  of  finding  the  A-T  factor  was  employed.  It  should  be  noted  that 
only  a  cursory  inspection  of  the  catchments  was  made  so  that  an  A-T  factor 
was  developed  for  the  study  area  as  a  whole  and  not  on  an  individual  water¬ 
shed  basis.  A  more  detailed  look  at  the  catchments  will  be  required  for  a 
more  accurate  definition  of  the  A-T  factor.  Table  A3  is  a  guide  for  selecting 
values  for  different  subfactors  used  to  obtain  the  A-T  factor.  Assigned  values 
for  each  category  are  summed  to  obtain  the  A-T  factor.  The  unit  debris  yield 
is  calculated  from  the  appropriate  equation  and  then  multiplied  by  the  A-T 
factor,  producing  an  adjusted  unit  debris  yield. 

The  watershed  characteristics  used  to  determine  the  subfactor  values  are 
presented  in  Table  A4.  Owing  to  the  cursory  nature  of  the  catchment  inspec¬ 
tion  for  the  sediment  impact  assessment,  two  A-T  factors  were  determined, 
one  A-T  factor  representing  an  average  case  and  one  a  case  with  more  severe 
erosion.  Assigned  subfactors  for  the  A-T  factor  determination  for  Fountain 
Creek  and  its  tributaries  can  be  seen  in  Table  A5.  The  tabulation  on  page  32 
gives  the  summary  results  of  sediment  yield  using  the  Los  Angeles  District 
method. 


Pacific  Southwest  Inter-agency  Committee  Method 
(PSIAC) 


This  technique  was  developed  to  evaluate  average  annual  sediment  yield  for 
a  variety  of  conditions  found  in  the  Pacific  southwest.  It  is  intended  for  broad 
planning  purposes  rather  than  specific  projects  where  more  intensive  investiga¬ 
tions  of  sediment  yield  would  be  required.  A  minimum  drainage  area  of 
10  square  miles  has  been  suggested,  but  for  the  reconnaissance  level  study  it 
has  been  applied  to  all  the  drainage  areas  irrespective  of  size.  For  the  smaller 
drainage  areas  this  would  increase  the  error  associated  with  the  calculation. 

The  PSIAC  method  places  the  area  under  consideration  in  one  of  five  sedi¬ 
ment  yield  classifications  using  nine  determining  factors.  The  nine  factors  are 
geology,  soils,  climate,  runoff,  topography,  ground  cover,  land  use,  upland 
erosion,  and  channel  erosion  and  sediment  transport.  Characteristics  of  each 


Chapter  3  Total  Sediment  Yield 


of  these  factors  can  be  seen  in  Table  A6,  with  an  appropriate  rating  alongside 
each.  Each  of  the  characteristics  has  been  assessed  for  its  contribution  to 
sediment  yield.  A  high  rating  is  given  to  a  characteristic  that  prompts  sedi¬ 
ment  yield,  and  so  on  for  a  moderate  and  low  rating.  Each  rating  has  been 
assigned  a  numerical  value  that  indicates  its  relative  significance.  The  sum  of 
these  numerical  values  for  the  nine  determining  factors  results  in  the  yield 
rating. 


Total  Sediment  Yield,  Los  Angeles  District  Method 

Stream 

Flood  Event 
(Percent 
Chance 
Exceedance) 

Drainage 

Area 

square 

miles 

Va 

Fire 

Factor 

riables 

A-T 

Factor 

Unit 

Debris 
-  Yield 
cu  yd/ 
square 

mile 

Unnamed 

1 ,  with 
no  burn 

1.4 

3 

1 

15,000 

21,000 

20,000 

28,000 

Black  Canyon  Creek 

1 ,  with 
no  burn 

2.6 

3 

0.55 

0.75 

9,600 

13,000 

25,000 

34,000 

Williams  Canyon  Creek 

1 ,  with 
no  burn 

2.7 

3 

0.55 

0.75 

9,900 

13,000 

27,000 

36,000 

Sutherland  Creek 

1,  with 
no  burn 

5.2 

3 

0.55 

0.75 

13,000 

17,000 

66,000 

90,000 

Ruxton  Creek 

1 ,  with 
no  burn 

12.1 

3 

10,000 

14,000 

120,000 

170,000 

Fountain  Creek 
(Equation  6) 

1 ,  with 
no  burn 

50.3 

3 

0.55 

0.75 

4,000 

5,000 

200,000 

270,000 

Fountain  Creek 
(Equation  7) 

1 ,  with 
no  burn 

50.3 

3 

0.55 

0.75 

3,300 

4,600 

170,000 

230,000 

Unnamed 

50,  with 

100%  burn 

1 .4 

6.5 

0.55 

0.75 

21,000 

28,000 

28,000 

38,000 

Black  Canyon  Creek 

50,  with 

100%  Burn 

2.6 

6.5 

0.55 

0.75 

13,000 

18,000 

34,000 

46,000 

Williams  Canyon  Creek 

50,  with 

1  00%  burn 

2.7 

0.55 

0.75 

13,000 

18,000 

36,000 

49,000 

Sutherland  Creek 

50,  with 

1  00%  burn 

5.2 

6 

0.55 

0.75 

1,400 

2,000 

7,500 

10,000 

Ruxton  Creek 

50,  with 

1  00%  burn 

12.1 

6 

0.55 

0.75 

1,100 

1,400 

13,000 

17,000 

Fountain  Creek 
(Equation  6) 

50,  with 

1  00%  burn 

50.3 

6 

0.55 

0.75 

370 

510 

19,000 

25,000 

Fountain  Creek 
(Equation  7) 

50,  with 

100%  burn 

50.3 

6 

0.55 

0.75 

170 

240 

8,700 

12,000 
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A  check  procedure  is  incorporated  into  the  technique.  A  high  sum  of  values 
A  through  G,  in  Table  A6,  should  result  in  high  summed  values  of  H  and  I. 
If  this  is  not  the  case,  then  either  unusual  erosion  conditions  exist  or  factors 
A-G  should  be  reevaluated.  The  conversion  of  the  yield  rating  is  done 
through  a  table  that  links  the  rating  to  a  class  and  an  average  annual  yield  as 
shown  in  the  following  tabulation. 


Rating 

Classification 

Sediment  Yield 
acre-feet/square  mile 

>100 

1 

>3.0 

75-100 

2 

1.0  -  3.0 

50  -  75 

3 

0.5  -  1.0 

25  -  50 

4 

0.2  -  0.5 

0  -  25 

5 

<  0.2 

The  ratings  assigned  to  each  of  the  nine  determining  factors  for  this  sedi¬ 
ment  impact  assessment  can  be  seen  in  Table  A7.  This  technique  was  carried 
out  with  only  a  cursory  viewing  of  the  watersheds,  so  that  although  individual 
subbasin  evaluation  was  possible,  it  was  done  on  an  elementary  level.  For  a 
more  detailed  and  accurate  employment  of  this  technique  a  more  detailed 
investigation  of  the  catchments  would  be  required.  More  error  will  be  in¬ 
herent  in  the  calculations  for  the  drainage  areas  that  are  less  than  10  square 
miles  as  they  are  smaller  than  the  limit  imposed  by  this  technique.  The  results 
for  the  average  annual  total  sediment  yield  follow: 


Stream 

Rating 

Total  Sediment  Yield 

PSIAC  Method 
acre-feet/square  mile 

Unnamed  Creek 

45 

0.2  -  0.5 

Black  Canyon  Creek 

35 

0.2  -  0.5 

Williams  Canyon  Creek 

42.5 

0.2  -  0.5 

Sutherland  Creek 

45 

0.2  -  0.5 

Ruxton  Creek 

45 

0.2  -  0.5 

Fountain  Creek 

37.5 

0.2  -  0.5 

Reservoir  Survey  Comparisons 

Reservoir  surveys  in  neighboring  watersheds  can  be  used  to  obtain  general¬ 
ized  regional  average  annual  sediment  yields.  Results  are  gathered  by  the 
Federal  Interagency  Subcommittee  on  Sedimentation  in  several  publications 
entitled  “Sediment  Deposition  in  U.S.  Reservoirs”  (U.S.  Department  of 
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Agriculture  1978;  Subcommittee  on  Sedimentation  1983,  1992).  Records  for 
central  and  southern  Colorado  and  the  northern  part  of  New  Mexico  between 
1960  and  1985  provide  some  basis  for  regional  comparisons.  Watersheds 
associated  with  the  reservoir  surveys  are  not  necessarily  similar  to  Fountain 
Creek  watershed.  Therefore,  the  regional  sediment  yields  should  be  used  only 
as  a  guideline.  Several  small  catchments  between  0.18  and  0.44  square  mile 
were  measured  on  Fort  Carson  Military  Base,  Colorado.  The  period  of 
record  varied  between  25  and  38  years.  The  average  annual  sediment  yield 
range  for  these  reservoirs  over  the  period  of  record  was  between  0.05  and 
0.57  acre-foot/square  mile.  Mud  Gulch  Reservoir  near  Canon  City,  CO,  has 
a  catchment  of  2.2  square  miles  and  has  an  annual  sediment  yield  of 
0.11  acre-foot/square  mile  over  a  10-year  period  of  record.  The  Abiquiu 
Reservoir  in  New  Mexico  is  large,  with  a  catchment  of  1,254  square  miles, 
and  has  an  average  annual  sediment  yield  over  a  15-year  period  of  record  of 
0.7  acre-foot/square  mile. 

Average  annual  sediment  yields  for  the  study  catchments  determined  using 
the  PSIAC  method  were  0.2  to  0.5  acre-foot/square  mile.  These  values  are 
within  the  range  of  data  from  the  regional  reservoir  sediment  yield  surveys. 

Total  Sediment  Yield  Results 


The  sediment  volumes  predicted  by  the  Los  Angeles  District  and  PSIAC 
methods  are  compared  in  the  following  tabulation: 


Stream 

Calculated  Sediment  Yield 

LA  District  Method 

PSIAC  Method 

Average  Annual 

Sediment  Volume 
cu  yd 

Sediment  Volume 

1  percent  Exceedance 

IMo  Burn 
cu  yd 

Sediment  Volume 

50  percent 

Exceedance 

100%  Burn 
cu  yd 

A-T 

Factor 

0.55 

A-T 

Factor 

0.75 

A-T 

Factor 

0.55 

A-T 

Factor 

0.75 

Unnamed 

20,000 

28,000 

28,000 

38,000 

450-  1,100 

Black  Canyon  Creek 

25,000 

34,000 

34,000 

46,000 

840  -  2,100 

Williams  Canyon  Creek 

27,000 

49,000 

870  -  2,200 

Sutherland  Creek 

66,000 

90,000 

7,500 

Hfljlifl 

1,700  -  4,200 

Ruxton  Creek 

120,000 

170,000 

17,000 

3,900  -  9,800 

Fountain  Creek 
(Equation  6) 

200,000 

270,000 

25,000 

16,000  -  41,000 

Fountain  Creek 
(Equation  7) 

170,000 

230,000 

12,000 

34 
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Sediment  yields  calculated  using  the  Los  Angeles  District  method  with  the 
lower  A-T  factor  can  be  considered  as  an  average  estimate.  The  yields  cal¬ 
culated  using  the  higher  A-T  value  should  be  considered  as  an  upper  range 
solution.  A  greater  range  in  results  may  be  obtained  using  the  standard 
deviation  associated  with  each  regression  equation.  It  needs  to  be  repeated 
that  only  a  preliminary  assessment  of  the  A-T  factor  was  possible  for  the 
sediment  impact  assessment  and  that  the  resultant  volumes  need  to  be  viewed 
with  this  in  mind.  As  Fountain  Creek  is  on  the  boundary  of  the  drainage  area 
restraints  for  two  equations,  both  have  been  used.  This  results  in  a  low  and 
high  sediment  volume  for  each  A-T  factor,  or  an  average  may  be  considered. 
It  is  interesting  to  note  that  the  smaller  catchments,  Unnamed,  Black  Canyon, 
and  Williams  Canyon  Creeks,  have  a  higher  50  percent  exceedance  volume 
than  the  1  percent  exceedance  volume.  This  is  the  direct  result  of  the  effects 
that  the  100  percent  burn  has  on  the  availability  of  debris  and  sediment. 
Catchments  that  have  experienced  no  burn  have  good  vegetation  so  that  the 
surfaces  are  protected.  The  larger  catchments,  Sutherland,  Ruxton,  and 
Fountain  Creeks,  have  much  lower  unit  sediment  yields  for  the  50  percent 
exceedance  flood,  even  with  100  percent  burn,  due  to  the  relatively  lower 
average  rainfall  that  occurs  over  the  larger  drainage  areas  during  the  50 
percent  exceedance  flood.  The  increase  in  sediment  potential  from  the  burn  is 
overshadowed  by  the  small  50  percent  chance  exceedance  unit  discharge. 
These  larger  catchments  also  show  lower  predicted  unit  sediment  yields  for 
the  1  percent  chance  exceedance  flood.  This  agrees  with  other  studies  that 
show  catchments  with  larger  drainage  areas  generally  exhibit  a  smaller 
delivery  ratio  than  smaller  ones.  Larger  catchments  will  usually  have  a  lower 
overall  slope,  smaller  upland  sediment  sources,  and  more  opportunity  for 
deposition,  all  of  which  reduces  the  potential  sediment  yield. 

Most  of  the  sediment  volume  produced  from  the  Unnamed,  Williams 
Canyon,  and  Ruxton  Creeks  will  not  reach  Fountain  Creek  because  pipes  and 
culverts  carry  the  flow  of  these  creeks  upstream  from  their  confluence  with 
Fountain  Creek.  These  pipes  constrict  the  channels  and  will  force  flood  flow 
into  the  streets.  Bed  load  remaining  in  the  creeks  will  deposit,  further  reduc¬ 
ing  conveyance  in  the  creek  and  causing  additional  flooding.  Sutherland 
Creek  has  a  constrictive  culvert  just  upstream  from  its  confluence  with 
Fountain  Creek  that  would  stop  the  majority  of  the  sediment  coming  down 
Sutherland  Creek  and  send  the  flow  out  into  the  streets.  Black  Canyon  Creek 
is  the  only  tributary  that  could  deliver  sediment  unobstructed  into  Fountain 
Creek. 

A  more  detailed  watershed  investigation  is  required  to  give  a  more  reliable 
estimate  of  the  sediment  yield  calculated  using  both  the  Los  Angeles  District 
and  the  PSIAC  Methods  during  the  next  level  of  planning  study.  A  detailed 
field  survey  would  be  required  of  each  subbasin.  This  would  involve  several 
different  areas  where  expertise  would  be  required:  geology,  soils,  geomor¬ 
phology,  river  hydraulics,  vegetation,  and  mountain  hydrology. 
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4  Debris  Flow  Potential 


Determining  Factors 


Debris  flows  are  the  main  process  ’•esponsible  for  the  formation  of  alluvial 
fans/debris  cones  at  the  base  of  mountain  valleys.  These  features  are 
generally  composed  of  a  poorly  sorted  mixture.  Debris  flows  are  a  form  of 
mass  wasting  and  occur  relatively  infrequently.  The  amount  of  fluid  in  a 
debris  flow  may  be  only  20  percent  or  less.  The  complementing  amount  of 
solid  material  can  vary  between  25  and  70  or  80  percent. 

The  occurrence  of  debris  flows  is  governed  by  several  factors: 

a.  The  composition  of  the  material  in  question  (there  needs  to  be  sufficient 
unconsolidated  debris  and  a  high  clay  contem). 

b.  The  underlying  beds  and  the  dip  of  the  bedding  planes. 

c.  The  steepness  of  slopes. 

d.  The  hydrological  situation. 

Conclusions  to  date  reveal  that  debris  flows  occur  on  sediment  beds  in  moun¬ 
tain  canyons  when  the  slopes  are  steeper  than  15  degrees  (0.27)  and  stop 
movement  when  the  slopes  are  less  than  3  degrees  (0.05).  The  finer  particles 
within  the  debris  flows  may  continue  movement  down  to  flatter  slopes  as  bed 
load  and  suspended  load.  Progressively  smaller  basins  with  steep  slopes  are 
more  suitable  to  debris  flows  than  their  larger  counterparts  because  thunder¬ 
storms  deliver  proportionally  larger  volumes  of  water  on  smaller  basins  and 
their  steeper  side  slopes,  resulting  in  instability  of  surficial  deposits,  making 
them  ready  for  transport. 

Gagoshidze  (1969)  summarized  the  requirements  for  debris  flows  in  the 
Soviet  Union: 

a.  Basins  should  have  a  predominance  of  rocks  rich  in  clay  forming 

aluminosilicates  and  clay  minerals  (clay,  argillites,  clay  shale,  gneisses, 
granitoid  rocks,  volcanic  ash,  tuff,  polymictic  sandstone,  etc.). 
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b.  The  basins  should  be  small,  <0.62  square  mile  per  mile  of  the  length 
of  the  main  watercourse. 

c.  The  main  watercourse  should  be  short,  <  18.6  miles  to  the  debris  cone. 

d.  The  basin  slopes  should  be  about  40  degrees  or  more,  the  average  slope 
of  the  main  watercourse  to  the  deposit  site  should  be  at  least  0. 10,  the 
channel  slope  of  any  lateral  tributaries  should  be  greater  than  the  main 
watercourse,  and  the  first  part  of  the  debris  cone  should  have  slopes  of 
at  least  0.04-0.05. 

e.  There  should  be  sharp  delineation  between  the  source  areas,  transport 
reaches,  and  deposition  sites. 

Antecedent  conditions  have  been  found  to  be  important.  A  rainfall  of 
0.25  in.  per  hour  in  an  area  where  the  total  seasonal  antecedent  rainfall  was 
10  in.  is  the  threshold  for  soil  slip  and  debris  flows  in  the  Santa  Monica 
Mountains  of  Southern  California.  Costa  and  Jarrett  (1981)  concluded  that 
debris  flows  can  occur  during  most  rainfall  events  provided  sufficient  material 
of  a  poorly  sorted  nature  is  available  on  side  slopes  ready  for  movement. 

Availability  of  material  is  an  important  factor  affecting  the  frequency  of 
debris  flow  events.  If  large  amounts  of  material  have  recently  been  removed 
by  a  large  runoff  event,  then  the  basin  is  not  able  to  produce  debris  flows 
again  until  sufficient  time  has  passed  to  replenish  loose  material  for  transport. 
This  is  accomplished  through  the  natural  processes  of  weathering,  mass 
wasting,  and  landslides.  In  larger  basins  the  data  seem  to  suggest  that  debris 
flows  do  not  occur  as  frequently  near  the  downstream  end  of  a  basin.  Old 
debris  flow  deposits  may  be  future  sources  of  debris  flow  or  may  be 
remobilized  during  future  events. 

Another  factor  that  affects  the  potential  for  debris  flows  is  the  condition  of 
the  watershed.  If  a  burn  has  occurred  in  the  area  recently,  the  potential 
increases  as  there  is  an  increase  in  source  sites.  Extensive  logging,  cattle 
grazing,  and  road  construction  will  also  increase  the  possibility  for  debris 
flows.  Vegetation  cover  is  a  major  factor  in  the  prevention  of  debris  flows  by 
reducing  the  debris  accumulation  and  deceasing  the  peak  discharge.  Costa  and 
Jarrett  (1981)  suggested  that  the  timberline  elevation  of  7,000-8000  ft  in  the 
Rocky  Mountains,  Colorado,  may  be  the  divide  between  two  environments 
where  hydrologic  processes  may  differ.  Below  this  elevation,  intense  rainfalls 
over  large  areas  are  frequent  and  may  cause  large  waterfloods  with  sediment 
loads  that  are  quite  high.  Above  this  elevation,  intense  rainfalls  are  less  fre¬ 
quent  and  debris  flows  may  occur  as  well  as  waterfloods. 


Analysis  of  Debris  Flow  Potential 

Available  guidelines  describing  the  prerequisites  for  debris  flow  are  listed 
in  the  following  tabulation: 


Chapter  4  Debris  Flow  Potential 


Creek 

Rock  Type 

Size 

square 
miles  per 
mile 

Creek 

Length 

miles 

Creek 

Slope 

Tribu¬ 

tary 

Slope 

Steeper 

Debris 

Cone 

Slope 

Sharp 

Delineation 

Unnamed 

Creek 

Precambrian 

igneous 

and 

metamorphic 

rocks: 

granites,  schists, 
gneisses. 

East 

narrow 

zone 

of 

Paleozoic 

and 

Mesozoic 

sedimentary 

rocks. 

0.90 

1.50 

0.29 

None 

Yes 

Black 

Canyon 

Creek 

0.87 

2.98 

0.13 

Yes 

None 

Yes 

Williams 

Canyon 

Creek 

0.63 

4.25 

0.12 

Yes 

None 

Yes 

Sutherland 

Creek 

1.27 

4.06 

0.17 

Yes 

None 

Yes 

Ruxton 

Creek 

2.03 

5.93 

0.18 

Most 

Yes 

None 

Yes 

Fountain 

Creek 

3.92 

12.83 

0.05 

Most 

Yes 

None 

Yes 

Because  no  detailed  field  investigation  was  conducted  for  the  sediment  impact 
assessment,  several  assumptions  were  made  regarding  the  prerequisites  for 
debris  flow.  Firstly,  the  rock  type  for  the  entire  Fountain  Creek  drainage  area 
was  based  on  observations  at  a  few  sites  and  assumed  to  be  homogeneous 
throughout.  The  tributary  slopes  are,  on  the  whole,  steeper  than  the  water¬ 
course  slopes,  and  the  assumption  was  made  that  there  was  sharp  delineation 
between  the  source,  transport,  and  deposition  sites.  Few  debris  source  sites 
and  no  debris  cones  were  observed  during  the  cursory  field  investigation.  The 
channel  slope  and  that  of  the  tributaries  were  calculated  using  the  blue  lines, 
which  represents  a  definable  channel,  from  a  1:24,000  quad  map  and  the 
elevations  taken  at  the  highest  point  at  the  end  of  the  blue  line  and  the  lowest 
elevation  at  the  potential  debris  site. 

The  observed  rock  type  of  the  drainage  area  is  in  a  category  that  can  pro¬ 
duce  debris.  The  main  watercourses  of  the  tributaries  are  short  and  the  chan¬ 
nel  slopes  are  steep  enough  to  support  movement  of  debris  flows.  This  is  not 
the  case  for  Fountain  Creek  itself.  These  factors  suggest  that  debris  flows  are 
possible  in  the  tributary  watersheds;  however,  no  debris  cone  deposits  were 
observed  and  the  drainage  basins  are  larger  than  0.62  square  mile/mile.  The 
area  is  well  vegetated  and  no  burn  has  occurred  in  the  area.  These  factors 
suggest  debris  flows  are  possible  but  unlikely  under  present  conditions  unless 
a  severe  storm  event  occurs  or  the  area  undergoes  a  burn,  logging,  or  road 
construction.  A  detailed  field  reconnaissance  is  required  to  fully  assess  debris 
flow  potential. 
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5  Summary  and 
Recommendations 


A  sediment  impact  assessment  was  conducted  for  Fountain  Creek  through 
Manitou  Springs,  CO.  The  purpose  of  the  study  was  to  identify  the  mag¬ 
nitude  of  sediment  problems  that  might  be  associated  with  proposed  flood- 
control  projects  and  to  recommend  appropriate  sedimentation  studies  for  the 
next  level  of  planning  study.  This  was  accomplished  using  a  sediment  budget 
approach  to  evaluate  channel  stability  in  terms  of  aggradation  and  degradation 
and  a  total  sediment  yield  analysis  to  provide  sediment  volumes  for  sizing 
reservoirs  and/or  debris  basins.  The  sediment  impact  assessment  provides 
only  qualitative  results  appropriate  for  evaluating  alternatives  at  the 
reconnaissance  level  planning  study. 

Available  data  were  used  in  the  sediment  impact  assessment.  Hydraulics 
were  based  on  an  existing  HEC-2  backwater  model  prepared  by  the 
Albuquerque  District.  Reach-averaged  values  for  hydraulic  parameters  were 
obtained  using  the  SAM  hydraulic  design  package.  Hydrology  used  included 
the  1  percent  and  50  percent  chance  exceedance  floods,  as  determined  by  the 
Albuquerque  District,  and  a  flow  duration  curve  developed  from  a  stream 
gauge  on  Fountain  Creek.  Bed  material  gradations  were  based  on  two 
samples  taken  to  be  representative  of  the  surface  and  subsurface  of  Fountain 
Creek  in  the  study  reach. 

Bed  material  sediment  yield  was  calculated  using  the  flow  duration- 
sediment  discharge-rating  curve  method.  Bed  material  sediment  yield  is  the 
sediment  load  primarily  responsible  for  channel  stability.  Sediment  transport 
was  calculated  using  three  different  sediment  transport  equations  to  obtain  a 
reasonable  range  of  possible  sediment  transport  rates.  Calculated  sediment 
yields  at  the  upstream  end  of  the  study  reach  for  the  1  percent  chance 
exceedance  hydrograph  varied  between  4,500  and  10,200  cu  yd.  Calculated 
sediment  yields  for  the  50  percent  chance  exceedance  hydrograph  varied 
between  310  and  920  cu  yd.  Average  annual  yield  varied  between  90  and 
380  cu  yd.  These  values  are  very  approximate  due  to  the  limited  data 
defining  the  character  of  the  streambed  and  the  variability  of  hydraulic 
characteristics  along  the  creek. 
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The  sediment  budget  analysis  demonstrated  the  variability  in  sediment 
transport  capacity  through  Manitou  Springs.  A  supply  reach  at  the  upstream 
end  of  the  study  reach  was  assigned.  Fountain  Creek  was  assumed  to  be  in 
equilibrium  in  the  supply  reach,  so  that  sediment  load  was  equal  to  the  sedi¬ 
ment  transport  capacity .  Aggradation  and  degradation  tendencies  were  deter¬ 
mined  for  upper,  middle,  and  lower  reaches  of  Fountain  Creek.  Both  degra¬ 
dation  and  aggradation  were  predicted  through  Manitou  Springs  during  floods 
and  on  an  average  annual  basis.  Between  52  and  36  percent  of  the  sediment 
load  delivered  from  the  supply  reach  should  deposit  in  the  upper  reach,  which 
includes  the  main  commercial  district  of  Manitou  Springs.  This  deposition  is 
caused  by  loss  of  sediment  transport  capacity  due  to  backwater  from  constric¬ 
tions  and  channel  obstructions.  The  deposition  will  reduce  channel  con¬ 
veyance  and  cause  maintenance  problems.  In  the  middle  reach,  where  there 
are  fewer  constrictions  and  channel  obstructions,  the  sediment  budget  analysis 
predicted  degradation.  This  will  occur  because  sediment  is  trapped  in  the 
upper  reach,  reducing  sediment  inflow  to  the  middle  reach  to  a  point  where  it 
is  less  than  the  sediment  transport  capacity.  Finally,  in  the  lower  reach, 
aggradation  is  predicted  due  to  reduction  in  channel  slope.  The  sediment 
budget  analysis  provides  an  estimate  of  the  magnitude  of  aggradation  and 
degradation  problems.  These  data  should  be  included  when  determining 
maintenance  costs  for  both  the  existing  channel  and  proposed  alternatives. 

One  flood-control  alternative  was  evaluated  for  its  impact  on  sedimentation 
in  Fountain  Creek.  Diverting  all  flood  flows  over  500  cfs  just  downstream 
from  Park  Avenue  would  result  in  significantly  more  deposition  in  Fountain 
Creek  downstream  from  the  diversion  in  the  commercial  district  of  Manitou 
Springs.  In  addition,  much  greater  degradation  was  predicted  downstream 
from  U.S.  Highway  24  where  the  diverted  flow  would  be  returned  to  Fountain 
Creek.  This  does  not  necessarily  render  the  proposal  infeasible,  but  additional 
maintenance  costs  will  be  expected. 

The  sediment  impact  assessment  included  calculated  estimates  of  total  sed¬ 
iment  yield  for  each  watershed  within  the  study  area.  Total  sediment  yield 
includes  the  wash  load  supplied  from  the  watershed  and  the  bed  material  load, 
which  is  governed  by  the  characteristics  of  the  streambed  and  the  flow. 

Average  annual  total  sediment  yield  and  the  total  sediment  yield  from  two 
storm  events  were  calculated.  The  storm  events  included  the  1  percent  and 
the  50  percent  chance  exceedance  floods.  The  50  percent  chance  exceedance 
flood  was  assumed  to  take  place  over  a  completely  burned  watershed.  The 
calculated  sediment  volumes  may  be  used  to  size  reservoirs  and/or  debris 
basins  during  the  reconnaissance  level  planing  study.  Maximum  calculated 
total  sediment  yields  for  the  1  percent  and  the  50  percent  chance  exceedance 
floods  are  13,000  and  28,000  cu  yd/square  mile,  respectively.  The  highest 
yield  rate  was  for  the  smallest  drainage  area.  Maximum  calculated  total 
sediment  yield  for  the  storm  events  considered  from  Fountain  Creek  ranged 
between  4,600  and  5,000  cu  yd/square  mile.  Average  annual  sediment  yield 
was  calculated  to  range  between  0.2  and  0.5  acre-foot/square  mile  for 
Fountain  Creek  and  all  the  tributaries. 
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The  potential  for  debris  flow  from  the  tributaries  was  investigated  for  the 
sediment  impact  assessment.  Although  some  of  the  watershed  characteristics 
that  produce  debris  flows  were  identified  in  the  Fountain  Creek  tributaries, 
there  was  no  evidence  of  historical  debris  flows  observed  during  the  limited 
field  reconnaissance.  The  preliminary  conclusion  is  that  debris  flows  are 
possible,  but  unlikely  unless  the  watershed  is  burned  or  a  severe  rainfall  event 
occurs.  Future  planning  studies  should  include  a  more  detailed  study  by  a 
geologist  and  geomorphologist. 

Due  to  the  variability  in  the  sediment  transport  capacity  in  Fountain  Creek, 
even  the  existing  channel  will  be  unstable  during  flood  conditions.  A  more 
detailed  sediment  study  is  recommended  for  the  next  level  of  planning  study. 
The  detailed  study  should  include  a  thorough  channel  bed  inventory, 
identifying  both  surface  and  subsurface  gradations  and  bedrock  outcrops  along 
Fountain  Creek.  It  is  expected  that  the  bed  material  gradations  will  vary 
longitudinally  through  the  study  reach.  Special  care  should  be  taken  to 
identify  bed  conditions  at  locations  where  degradation  is  predicted.  Locations 
and  reliability  of  channel  bed  stabilizers  should  also  be  determined.  The 
cross-section  definition  in  the  HEC-2  backwater  model  used  in  this  study  is 
generally  inadequate  for  a  feasibility  level  channel  stability  or  sedimentation 
study.  More  detailed  surveys  of  cross-section  geometry  will  be  required  at  the 
next  level  of  planning  study.  New,  more  detailed  topographic  mapping  is 
available  and  will  provide  better  geometric  definition.  However,  field  surveys 
will  still  be  required  to  supplement  the  topographic  data.  Much  better 
quantitative  estimates  for  sediment  deposition  and  scour,  for  both  existing 
conditions  and  alternative  designs,  can  be  obtained  using  the  HEC-6  numerical 
sedimentation  model.  This  model  accounts  for  variable  armoring  and  for  the 
variation  in  channel  cross-sectional  areas.  Careful  adjustment  of  the  model 
would  produce  a  more  reliable  estimate  of  bed  material  transport  through 
Manitou  Springs.  The  model  should  be  adjusted  to  reasonably  replicate 
existing  conditions,  i.e.,  relatively  small  changes  during  a  50  percent  chance 
exceedance  hydrograph.  A  reliable  sediment  transport  equation  would  need  to 
be  identified.  Measured  sediment  transport  data  from  similar  streams  in  the 
region  could  be  used  to  select  a  sediment  transport  equation. 

The  sediment  yield  estimates  were  calculated  using  very  limited  watershed 
data.  If  a  storage  option  is  pursued  in  the  next  level  of  planning  study,  a 
much  more  detailed  investigation  of  the  watershed  will  be  required  and  more 
refined  calculations  made.  Input  for  the  more  detailed  sediment  yield  study 
should  be  obtained  from  several  disciplines  including  geology,  soils, 
hydraulics,  and  hydrology. 

The  following  tasks  are  identified  for  a  more  detailed  sedimentation  study 
of  Fountain  Creek: 

a.  Channelization/Diversion  Alternatives 
(1)  Channel  bed  inventory 
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(a)  Bed  material  sampling  program 

(b)  Location  of  hard  points  and  bedrock  outcrops 

(c)  Location,  dimensions,  and  stability  of  grade  control  structures 

(2)  Development  of  geometric  model 

(a)  Cross-section  layout  based  on  locations  of  grade  control 

(b)  Field  surveys 

(3)  Selection  of  sediment  transport  equation 

(a)  Locate  and  analyze  measured  sediment  data  from  similar 
stream 

(b)  Evaluate  sediment  transport  equations 

(4)  Develop  HEC-6  model  of  existing  conditions 

(a)  Adjust  model  to  simulate  existing  conditions  for  50  percent 
chance  exceedance  flood 

(b)  Simulate  1  percent  chance  exceedance  flood 

(c)  Long-term  simulation 

(5)  Evaluate  alternatives 

(a)  Simulate  50  and  1  percent  chance  exceedance  floods 

(b)  Perform  a  long-term  simulation  to  determine  maintenance 
requirements 

b.  Storage  Alternatives 

(1)  Sediment  yield  from  watershed 

(a)  Determine  watershed  characteristics 

(b)  Use  at  least  two  methods  to  calculate  average  annual  sediment 
yields 

(2)  More  detailed  study  for  debris  flow  potential 


Chapter  5  Summary  and  Recommendations 
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Appendix  A  Sediment  Yield  Tables 


Table  A2 

Calculated  Sediment  Yields,  All  Reaches 


_ -II 

Average  Annual 

50  Percent  Chance 

Exceedance 

1  Percent  Chance 
Exceedance 

cubic 

yards 

Percent  of 
Upstream 

cubic 

yards 

Percent  of 
Upstream 

cubic 

yards 

Percent  of 
Upstream 

Supply  Reach 

Yang 

150 

660 

10,200 

Ackers- 

White 

93 

310 

4,500 

Toffaleti- 

Schoklitsch 

380 

920 

7,100 

Reach  2 

Yang 

85 

57 

420 

64 

6,400 

63 

Ackers- 

White 

42 

45 

190 

61 

2,900 

64 

Toffaleti- 

Schoklitsch 

157 

41 

580 

62 

4,700 

66 

Reach  3 

Yang 

93 

109 

410 

98 

7,600 

118 

Ackers- 

White 

44 

105 

180 

97 

3,400 

116 

Toffaleti- 

Schoklitsch 

150 

96 

570 

98 

5,900 

126 

Reach  4 

Yang 

150 

160 

680 

166 

20,000 

260 

Ackers- 

White 

81 

184 

310 

170 

8,600 

255 

Toffaleti- 

Schoklitsch 

320 

213 

990 

174 

18,000 

305 

Reach  5 

Yang 

110 

74 

560 

82 

9,500 

48 

Ackers- 

White 

55 

68 

250 

80 

4,200 

48 

Toffaleti- 

Schoklitsch 

200 

63 

790 

80 

9,600 

53 
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Table  A3 

Los  Angeles  District  A-T  Factor  Table 

A-T  Subfactor 

0.25 

0.20 

0.15 

0.10 

0.05 

Subfactor  Group  1:  Parent  Material 

Folding 

Severe 

Moderat 

e  to 

Severe 

Moderate 

Minor  to 

Moderate 

Minor 

Faulting 

Severe 

Moderate 

Minor 

Fracturing 

Severe 

Moderate 

Minor 

Weathering 

Severe 

Moderate 

Minor 

Subfactor  Group  2:  Soils 

Soils 

Noncohesive 

Partly 

Cohesive 

Highly 

Cohesive 

Soil  Profile 

Minimal  Soil 
Profile 

Some  Soil 
Profile 

Well- 

Developed 

Soil  Profile 

Soil  Cover 

Much  Bare 

Soil  in 

Evidence 

Some  Bare 
Soil  in 
Evidence 

Little  Bare 

Soil  in 

Evidence 

Clay 

Colloids 

Few  Clay 
Colloids 

Some  Clay 
Colloids 

Many  Clay 
Colloids 

Subfactor  Group  3:  Channel  Morphology 

Bedrock 

Exposures 

Few 

Segments  in 
Bedrock 

Some 
Segments 
in  Bedrock 

Many 

Segments 
in  Bedrock 

Bank 

Erosion 

>30%  of 
Banks 

Eroding 

10-30%  of 
Banks 

Eroding 

<10%  of 

Banks 

Eroding 

Bed  and 

Bank 

Material 

Noncohesive 
Bed  and 

Banks 

Partly 

Cohesive 

Bed  and 

Banks 

Highly 

Cohesive 

Bed  and 

Banks 

Vegetation 

Poorly 

Vegetated 

Some 

Vegetation 

Much 

Vegetation 

Headcutting 

Many 

Headcuts 

Few 

Headcuts 

No 

Headcuts 

(Continued) 

Sum  all  A-T  Subfactors  from  the  four  groups  =  A-T  Factor 
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Table  A3  (Concluded) 

A-T  Subfactor 

0.25 

0.20 

0.15 

0.10 

0.05 

Subfactor  Group  4:  Hill  Slope  Morphology 

Rills  and 

Gullies 

Many  and 
Active 

Some  Signs 

Few  Signs 

Mass 

Movement 

Many  Scars 
Evident 

Few  Signs 
Evident 

No  Signs 
Evident 

Debris 

Deposits 

Many 

Eroding 

Deposits 

Some 

Eroding 

Deposits 
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Table  A4 

Input  Parameters  for  Los  Angeles  District  Method 
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Table  A5 

A-T  Factor  Determination 

A-T  Subfactor 

Average  Case 

More  Severe  Case 

Parent  material 

0.15 

0.15 

Soils 

0.15 

0.20 

Channel  morphology 

0.15 

0.20 

Hill  slope  morphology 

0.10 

0.20 

A-T  factor 

0.55 

0.75 
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Table  A6 

Factors  Affecting  Sediment  Yield,  PSIAC  Method 

Sediment 

Yield 

Levels 

A 

Surface 

Geology 

B 

Soils 

C 

Climate 

D 

Runoff 

E 

Topography 

High 

(10) 

(10) 

(10) 

(10) 

(20) 

a.  Marine 

shales  and 

related 

mudstones  and 

siltstones. 

a.  Fine  textured; 
easily  dispersed; 
saline-alkaline; 
high  shrink- 
swell 

characteristics. 

b.  Single  grain 
silts  and  fine 

sands 

a.  Storms  of 
several  days' 
duration  with 
short  periods 
of  intense 

rainfall 

b.  Frequent 
intense 

convective 

storms 

c.  Freeze-thaw 

occurrence 

a.  High  peak 
flows  per  unit 

area 

b.  Large 
volume  of 
flow  per  unit 

area 

a.  Steep  upland 
slopes  (in  excess 
of  30%);  high 
relief;  little  or  no 
floodplain 
development 

*  * 

Moderate 

(5) 

(5) 

(5) 

(5) 

(10) 

a.  Rocks  of 

medium 

hardness 

b.  Moderately 
weathered 

c.  Moderately 
fractured 

a.  Medium 

textured  soil 

b.  Occasional 
rock  fragments 

c.  Caliche  layers 

a.  Storms  of 

moderate 
duration  and 
intensity 

b.  Infrequent 
convective 

storms 

a.  Moderate 
peak  flows 

b.  Moderate 

volume  of 
flow  per  unit 

area 

a.  Moderate 
upland  slopes 
(less  than  20%) 

b.  Moderate  fan 
or  floodplain 
development 

*  * 

Low 

(0) 

a.  Massive, 
hard 

formations 

(0) 

a.  High 
percentage  of 
rock  fragments 

b.  Aggregated 
clays 

c.  High  in 
organic  matter 

(0) 

a.  Humid 
climate  with 
rainfall  of  low 
intensity 

b.  Precipitation 
in  form  of 

snow 

c.  Arid  climate, 
low-  intensity 
storms 

d.  Arid  climate; 
rare  convective 

storms 

(0) 

a.  Low  peak 
flows  per  unit 

area 

b.  Low 

volume  of 
runoff  per 
unit  area 

c.  Rare  runoff 

events 

(0) 

a.  Gentle  upland 

!  slopes  (less  than 

5%) 

b.  Extensive  allu¬ 
vial  plains 

(Continued) 

Note  -  The  numbers  in  specific  boxes  are  the  values  to  be  assigned  to  the  characteristics, 

**  If  experience  so  indicates,  interpolation  between  the  three  sediment  yield  levels  may  be  made. 

A8 


Appendix  A  Sediment  Yield  Tables 


Table  A6  (Concluded) 


Sediment  Yield 

Levels 

F 

Ground 

Cover 

G 

Land  Use 

H 

Upland 

Erosion 

1 

Channel  Erosion 
and  Sediment 

Transport 

High 

{10} 

Ground  cover 
does  not  exceed 
20% 

a.  Vegetation 
sparse;  little  or 
no  litter 

b.  No  rock  in 
surface  soil 

(10) 

a.  More  than  50% 
cultivated 

b.  Almost  all  of 
area  intensively 
grazed 

c.  All  of  area 
recently  burned 

(25) 

a.  More  than  50% 
of  the  area 
characterized  by  rill 
and  gully  or 
landslide  erosion 

(25) 

a.  Eroding  banks  con¬ 
tinuously  or  at 
frequent  intervals  with 
large  depths  and  long 
flow  duration 

b.  Active  headcuts 
and  degradation  in 
tributary  channels 

*  * 

Moderate 

(0) 

Cover  not 
exceeding  40% 

a.  Noticeable 
litter 

b.  If  trees 
present,  under¬ 
story  not  well 
developed 

(0) 

a.  Less  than  25% 
cultivated 

b.  50%  or  less 
recently  logged. 

c.  Less  than  50% 
intensively  grazed 

d.  Ordinary  road 
and  other 
construction 

(10) 

a.  About  25%  of 
the  area 

characterized  by  rill 
and  gully  or 
landslide  erosion 

b.  Wind  erosion 
with  deposition  in 
stream  channels 

(10) 

a.  Moderate  flow 
depths,  medium  flow 
duration  with 
occasionally  eroding 
banks  or  bed 

*  # 

Low 

(-10) 
a.  Area 

completely  pro¬ 
tected  by 
vegetation,  rock 
fragments,  litter. 
Little 

opportunity  for 
rainfall  to  reach 
erodible  materia! 

(-10) 

a.  No  cultivation 

b.  No  recent 
logging 

c.  Low-intensity 
grazing 

(0) 

a.  No  apparent 
signs  of  erosion 

(0) 

a.  Wide  shallow 
channels  with  flat 
gradients,  short  flow 
duration 

b.  Channels  in  massive 
rocks,  large  boulders 
or  well  vegetated 

c.  Artificially 
controlled  channels 
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Table  A7 

Assigned  Ratings,  PSIAC  Method 

Determining 

Factor 

Un¬ 

named 

Creek 

Black 

Canyon 

Creek 

Williams 

Canyon 

Creek 

Sutherland 

Creek 

Ruxton 

Creek 

Fountain 

Creek 

A.  Surface 

Geology 

5 

5 

5 

5 

5 

5 

B.  Soils 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

C.  Climate 

7.5 

7.5 

7.5 

7.5 

7.5 

7.5 

D.  Runoff 

7.5 

7.5 

7.5 

7.5 

7.5 

7.5 

E.  Topography 

20 

10 

17.5 

20 

20 

12.5 

F.  Ground  Cover 

*7.5 

-7.5 

-7.5 

-7.5 

-7.5 

-7.5 

G.  Land  Use 

-10 

-10 

-10 

-10 

-10 

-10 

H.  Upland 

Erosion 

5 

5 

5 

5 

5 

5 

I.  Channel 
Morphology  and 
Sediment 

Transport 

15 

15 

15 

15 

15 

15 

Total  Rating 

45.5 

35 

42.5 

45 

45 

37.5 
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